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We report the results of theoretical investigation of electronic structure of ThF^ cation which is 
one of the most interesting systems to search for the permanent electric dipole moment (EDM) of 
the electron (eEDM) [H. Loh, K.C. Cossel, M.C. Grau, K.-K. Ni, E.R. Meyer, J.L. Bohn, J. Ye, 

E.A. Cornell, Science 342, 1220 (2013)] and other effects of violation of time reversal (T) and 
spacial parity (P) symmetries in fundamental interactions. For the working ®Ai state we have 
found a quite high value of the effective electric field acting on unpaired electrons (37.3 GV/cm). 

The field will be required to interpret the experiment planed on ThF"*" in terms of eEDM. Within the 
concept of atoms in compounds [A.V. Titov, Y.V. Lomachuk, and L.V. Skripnikov, Phys. Rev. A 
90, 052522 (2014)] we have compared the ThF+ electronic structure with that of ThO. Also we 
have calculated other parameters of T,P-odd interactions: Wt,p, which is needed for interpretation 
of the experiment in terms of the dimensionless constant kr.p characterizing the strength of the 
T,P-odd pseudoscalar—scalar electron—nucleus neutral current interaction (50 kHz); Wm, which is 
required to search for the Th nuclear magnetic quadrupole moment in ^^^ThF^ (0.88 ,.,^2 )■ A 

number of properties which can be measured are also calculated: hyperfine structure constant, the 
molecule-frame dipole moment, and g-factor. 

INTRODUCTION 

During the last decade an impressive progress in the 
search for the permanent electric dipole moment of the 
electron (eEDM) has been achieved [TH3]. The great in¬ 
terest in eEDM is caused by the fact that its nonzero 
value implies manifestation of interactions which are not 
symmetric with respect to both time (T) and space 
(P) inversions (T,P-odd interactions). According to the 
Standard model eEDM should be less than 10“^®e-cm [1]. 

Therefore the observation of eEDM at a notably higher 
level would indicate the presence of a “new physics” be¬ 
yond the Standard model. Most popular extensions of 
the Standard model predict the magnitude of the eEDM 
at the level of 10“^® — 10“^®e-cm m and that range is 
almost passed to-date by the latest atomic and molecular 
measurements. 

It was found since sixties of the past century HIMl] 
that extremely sensitive experiments towards the search 
of T,P-odd effects can be performed on heavy-atom 
molecules and solids. The current limit, \de\ < 8.7x10“^® 
e-cm (90% confidence), was set with a molecular beam 
of thorium monoxide (ThO) molecules in the metastable 
electronic state [5]. The previous best limit was 

also established on a molecular beam but using the YbF 
radicals [2]. 

Nowadays, a number of new prospective systems are 
suggested, investigated theoretically and, in part, pre¬ 
pared experimentally (HfF+ [T^HTT] . YbF ITBifel] . 

ThO [1 HSISS], ThF+ [g, WC [El IE], PbF lEl- 
155] .RaO [5S1 177] . RaF [El I3S] etc.) which promise to 
achieve a sensitivity to eEDM up to 10“^® — 10“®®e-cm. 


One of promising experiments towards the measurement 
of eEDM is proposed on the ®Ai state of the cation of 
thorium monofluoride (ThF+) by E. Cornell group [14] . 
The use of the ^ Ai state has a number of advantages from 
experimental point of view. Due to D-doublet structure 
of ® Ai state the interval between the opposite parity lev¬ 
els is very small. Therefore, the molecule can be polar¬ 
ized by a weak electric field which leads to cancellation 
of some systematic errors since the effect on the dou¬ 
blet components has an opposite sign [25l uni Ell . Also, 
magnetic moment (g-factor) of the ®Ai electron state is 
very small (zero in the nonrelativistic limit), and this is 
another reason for reducing the systematic errors. The 
advantage of using such a state has been demonstrated 
in the recent experiment on ThO molecule [3]. 

The working ^Ai state of ThO is a metastable (first 
excited) one with the lifetime of about 2 ms [E] whereas 
the ground state is In contrast to the ThO case, the 
energies of the ®Ai and states in ThF+ are very close 
[Ej . In Ref. ED the ®Ai state of ThF+ was assigned as 
the first excited state with transition energy 316 cm“^ 
(versus 5321 cm“^ in ThO [E]). However, the most re¬ 
cent experiments by Cornell group show that ®Ai is the 
ground state of ThF+ [55]. This suggests a very good 
statistics. 

To interpret the results of the ThF'*' experiment in 
terms of the eEDM one should know a parameter usu¬ 
ally called “the effective electric field on electron”. Egg, 
which cannot be measured. Egg actually relevant to only 
the spin-polarized electrons (the closed shells do not con¬ 
tribute to measured effects in context of eEDM, see next 
section), it can be evaluated as an expectation value of 
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the T,P-odd operator (see Refs. mHni): 




Hd{i) 


1 ^), 


( 1 ) 


where de is the value of eEDM,'!' is the wave function of 
the considered state of ThF+, and n = (d>|J • n|d>), J 
is the total electronic momentum, n is the unit vector 
along the molecular axis C directed from Th to F (fl = 1 
for the considered ^Ai state of ThF+), 

Hd = 2de f n E 1 . (2) 


^0 (tE J ’ 

E is the inner molecular electric field, and cr are the Pauli 
matrices. In these designations Egg = IFd|n|. 

Besides the interaction given by operator ([^ there 
is a T,P-odd pseudoscalar—scalar electron—nucleus neu¬ 
tral currents interaction with the dimensionless constant 
kT,p- Note that it was estimated in Ref. [SO] within the 
Standard model that this interaction can induce even 
greater T,P-odd effect in ThO simulating the eEDM. The 
interaction is given by the following operator [51) : 

Gp 

Ht,P = i^^ZkT,P^0j5PN{Tc), (3) 

where Gf is the Fermi constant, 70 and 75 are the Dirac 
matrixes and pn(j) is the nuclear density normalized to 
unity. To extract the fundamental kp^p constant from 
an experiment one needs to know an electronic structure 
factor, Wt,p, on a nucleus of interest: 


i. ’ 


(4) 


Similarly to E^g, the Wt,p parameter cannot be mea¬ 
sured and have to be obtained from a molecular electronic 
structure calculation. 

In Refs. [23 IS2] is was demonstrated that ^^®ThO 
molecule can be used to search for T,P-odd interaction 
of ^^®Th nuclear magnetic quadrupole moment (MQM) 
with electrons. The T,P-odd electromagnetic interaction 
is described by the Hamiltonian iniiTI |53|; 

=-(5) 

2/(2/- I) 2 A ’ ^ 

where Einstein’s summation convention is implied, a are 
the 4x4 Dirac matrices, 0 ;=^^ ^^,ris the displace¬ 
ment of the electron from the Th nucleus, I is the nuclear 
spin, M is the nuclear MQM, 

~ 2 /( 2 / - 

Ti,k = lilk + Ikli ~ |di,fed(/ + 1) . (7) 

In the subspace of ±D states Hamiltonian ([^ is reduced 
to the following effective molecular Hamiltonian |9]: 
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where S is the effective electron spin |^, 5=1/21=1. Wm 
parameter can be evaluated by the following matrix ele¬ 
ment [52] : 

It was shown in Refs. |2H1[S1] that using the ^^®Th isotope 
one can obtain limits on the strength constants of T,P- 
odd nuclear forces, neutron EDM, QCD vacuum angle 0, 
quark EDM and chromo-EDM. This also can be applied 
to ^^®ThF+ cation. 

A commonly used way of verification the theoretical 
Eeff, Wt,p and Wm values is to calculate “on equal foot¬ 
ing” (using the same approximation for the wave func¬ 
tion) those molecular characteristics (properties or ef¬ 
fective Hamiltonian parameters) which have comparable 
sensitivity to different variations of wave function but, in 
contrast, can be measured. Similar to Ef.g, Wt,p and 
Wm these parameters should be sensitive to a change in 
the spin-polarized share of the electronic density, etc., in 
the atomic core region. The hyperfine structure (HFS) 
constant, A||, is traditionally used as such a parame¬ 
ter (e.g., see Ref. jSS]). To obtain A|| on ^^®Th in the 
229ThF+ theoretically, one can evaluate the following ma¬ 
trix element: 



where /rxh is magnetic moment of an isotope of ^^®Th 
nucleus having spin I. In the present paper we do not 
consider fluorine nuclear spin. 

For the preparation and conduction of the experiment 
the value of g-factor of the molecule is of interest. It is 
defined as 

Gi\ = ^{^\Ll-gsSm, ( 11 ) 

where L® and 5® are the electronic orbital and electronic 
spin momenta operators, respectively; gs = —2.0023 is 
a free—electron ^-factor. Note that the value of Gy is 
close to zero for the ^Ai state (and equal to zero when 
both the scalar-relativistic approximation is applied and 
the radiation corrections to the free-electron g-factor are 
ignored. Therefore, the parameter is very sensitive to 
the quality of the wave function, since high-order inter¬ 
ference contributions between the spin-orbit and electron 
correlation effects become important. 

Recently, ThF+ has been studied both experimentally 
and theoretically in Refs. |43l |44]. The measured and 
calculated values are given there for spectroscopic con¬ 
stants of the lowest-lying states including ^Ai. The lat¬ 
ter is found to be the first excited state. However, up 
to now there is only one (semiempirical) estimate of E^g 
in ThF+ published in Ref. [26], Eeg=90 GV/cm. The 
aim of the paper is to perform accurate ab-initio study 
of ThF+ electronic structure and calculate E^g and other 
parameters given by Eqs. 
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TWO-STEP APPROACH 


It follows from Eqs. ([^-(|l0|) that the action of oper¬ 
ators related to the Eeff, Wt,p, Wm and A|| character¬ 
istics is heavily concentrated in the atomic core region. 
On the other hand the leading contribution to the cor¬ 
responding matrix elements (mean values) is due to the 
valence electrons since contributions from the inert (usu¬ 
ally closed and spherically symmetric) inner-core shells 
compensate each other or negligible in most cases of prac¬ 
tical interest for the operators, in particular, dependent 
on the total angular momentum and spin. Note, how¬ 
ever, that the spin-polarization of core (sub-valence or 
outer-core) shells induced by the valence unpaired elec¬ 
trons can provide a comparable contribution by magni¬ 
tude to such properties as that from the valence elec¬ 
trons, e.g., see Refs. [BISIIISS]- Below we shall call such 
properties as the “core properties” (or, more generally, 
“core characteristics” since not only measurable prop¬ 
erties but other effective Hamiltonian parameters which 
are not always measurable can be considered here) as¬ 
suming that the main contribution to them comes from 
the spatially-localized core region rather than from core 
shells. Some of well-known examples of such property 


are the magnetic dipole hyperfine constants (see Eq. 10). 
In the cases of unpaired s-electrons, the leading contri¬ 
bution to the hyperfine structure is determined by the 
Fermi-contact interaction (in nonrelativistic case) which 
is proportional to the electronic spin density directly on 
the nucleus. Other examples of core properties (which 
have negligible contribution from inert inner-most core 
shells) are the chemical shifts of X-ray emission spectra 
[S711SH], etc. 

One can safely exclude inactive inner-core electrons 
from correlation calculation due to their negligible con¬ 
tribution to the core properties. In the present consider¬ 
ation, the inner-core consists of Is—4/ electrons of Th. 
As Th is a very heavy element (atomic number is 90) 
the interaction of electrons with the Th nucleus should 
be treated by a fully relativistic manner for a good accu¬ 
racy. Moreover, for some properties even taking account 
of Breit interaction (mainly between valence and core 
electrons of Th, see [sniiso]) can be important. With 
a good accuracy for the properties considered here, the 
inner-core electrons differ negligibly in the cases of atomic 
Th and ThF+ cation because their wavefunctions mostly 
defined by the strong Th nucleus potential ~ ^ screened 
by inner-more electrons, so that the effective Th core field 
is much stronger than the energetics of valence (chemi¬ 
cally active) electrons. In the correlation calculation they 
can be frozen without a loss of the accuracy accessible 
presently. A common way to exclude inner-core electrons 
is to use the relativistic effective core potential method. 
Earlier our group has developed the generalized relativis¬ 
tic effective core potential (GRECP) version which per¬ 


mits one to attain a very high accuracy [STH55] . This ef¬ 
fective potential emulates interaction between inner-core 
electrons (excluded explicitly from GRECP calculations) 
and valence plus outer-core electrons (treated explicitly 
with GRECP). 

Performing electronic structure calculation one can 
evaluate different valence properties such as transition 
energies between low-lying states, molecule-frame dipole 
moments, etc. However, since the inner-core parts 
of the valence one-electron “pseudo-wavefunctions” are 
smoothed in the GRECP calculations, they have to be 
recovered with some core-restoration method before us¬ 
ing them to evaluate the core characteristics considered 
above. In series of papers a non-variational restoration 
concept (and its initial implementation, see [49] and ref¬ 
erences), which is based on a proportionality of valence 
and low-lying virtual spinors in the inner-core region of 
heavy atoms was developed (see [58] and the next sec¬ 
tion) . Recently we have developed a new implementation 
of the concept which permits to use well-developed codes 
on correlation treatment such as DIRAC (Mj, mrcc [S5] 
and CFOUR [nni [53 157] . Below we give description of 
the new implementation; now the code is also extended 
to characterize effective states (configurations) of atoms 
in compounds, e.g. Th in the ThF'*' and ThO. 

Using the basic idea of the nonvariational restoration 
method one generates equivalent basis sets of one-center 
four-component spinors 

f fnlji'f')9ijm \ 

\ 9nljir)02j-l,jm J 

and smoothed two-component pseudospinors 
fnlj (l ’^^Ijm 

in all-electron finite-difference Dirac-Fock-Breit and 
GRECP / self-consistent field calculations (employing the 
jj—coupling scheme) of the same configurations of a con¬ 
sidered atom and its ions EHHZI] [za. Here n is the prin¬ 
cipal quantum number, j is the total electronic momen¬ 
tum, m is its projection and I is the orbital momentum. 
In the newly developed procedure a basis set of real spin- 
orbitals (and not complex spin-orbit-mixed spinors) 
is generated additionally. The spin-orbitals are then 
expanded in the basis set of one-center two-component 
atomic pseudospinors 

i„^axi = l^ + l/2| 

~ ^ ' ( 1 ^) 

1=0 i=\l-\l2\n,m 


The atomic two-component pseudospinors are replaced 
by equivalent four-component spinors while the expan¬ 
sion coefficients from Eq. (12) are preserved: 




imaxi = |i + l/2| 

E E Ej’»i 


1=0 j = |;-l/2| n,m 


gnlj{'r)02j-ljm J 


(13) 
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and we obtain four-component function which is 
“equivalent” to ^p. If a one-electron reduced den¬ 
sity matrix with elements in a basis set of multi¬ 
center spinors (or spin-orbitals) ipp is evaluated after the 
(G)RECP calculation of a molecule or some condensed- 
matter system (see Ref. m for details on the condensed- 
matter case) one can then reexpand it in the basis of one- 
center ^p functions on an atom of interest. This mapping 
from a multi-center basis {ipfj.} to the one-center basis 
{^p} corresponds to a similarity transformation of the 
density matrix: 

||.P/il/|| ll^pqll (14) 


where Dpq are elements of the density matrix in the ba¬ 
sis of ^p functions. Due to “equivalence” of ^p and ^p 
functions (see Eqs. (12) and (Hi) based on appropriate 
properties of the hard-core shape-consistent (G)RECP 
versions |73] . one can write: 


for which the mean values of operators corresponding 
to all the considered core characteristics, have near the 
same magnitudes for the case of an atom bonded in a 
molecule and for the same atom in the considered ef¬ 
fective state. Note, that the concept cannot be directly 
applied to evaluation of the “valence” properties (again, 
taking in mind spatial localization rather than affiliation 
to valence shells) such as the molecule-frame dipole mo¬ 
ment, g-factor, etc. 

Let us show how the AIC theory can be formulated in 
context of the problems discussed in this paper. Assume 
that we have obtained a one-electron density matrix from 
calculation of an atom, molecule or crystal (in the direct 
lattice). The density matrix can be formally reexpanded 
on one center, i.e., on a heavy atom of interest: 

p(T|r') = 

nljm,n'l'j'm' 


\D. 


pg\ 


\D. 


pq\ 


(15) 


where Dpg are elements of the density matrix in the basis 


of four-component ^p functions (13). Thus, as an approx¬ 


imation we can equate the Dpq elements to Dpq and inter¬ 
pret it as a restoration of “true” four-component struc¬ 
ture of density matrix that is important first of all for 
the inner-core region. 

The mean value of some one-electron operator A cor¬ 
responding to a core property on a given atom can be 
evaluated as follows: 


(A) — ^ DpqAp 


(16) 


pq 


in a sufficiently complete basis set of orthonormal atomic 
functions {’Pnijm}- Then for the mean value of some one- 
electron operator A we have: 


{A) — Pnljm,n'l'j'm' J ‘T’l'/'j'm'' 

nljm,n'l' j'm' 

(18) 

The mean value in Eq. can be rewritten as 


I J ^1, p j'm' 

NI<|Ad 

|r|>|Ac| 


(19) 


where Apq are the matrix elements of operator A in the 
basis of four-component functions ^p (13). 

In the current implementation of restoration procedure 
the functions ^p are real spin-orbitals with the spatial fac¬ 
tor in the form of contracted Gaussians, for which reex¬ 
pansion (14) is performed analytically. Therefore, signif¬ 
icant acceleration is attained in contrast to the original 
restoration procedure [49]. Note, however, that the four- 
component functions used to evaluate matrix elements 
of operator A are taken in numerical (finite-difference) 
form. This permits to exclude some complications in 
reproducing accurate wavefunction behavior in a region 
near nucleus which can arise when Gauss-type functions 
are used there. 


“ATOMS IN COMPOUNDS” THEORY 

In Ref. |58j we have introduced a concept of atoms in 
compounds (AIC) and applied it to the problem of chem¬ 
ical shifts of X-ray emission lines. The concept assumes 
that for the core characteristics one can determine an ef¬ 
fective state of a given atom in a chemical compound. 


where Rc is some “core radius” (see also below). Here we 
consider that the operator A corresponds to a core prop¬ 
erty. It means that for r>Rc the second term in Eq. (19) 
have to be negligible compared to the first term. As an 
extremal case, Rc—0 for the Eermi-contact interactions. 
Thus, for a core-property operator we have 


\P\<\Ro\ 

( 20 ) 

Now we assume that the basis set {<Pnijm} was con¬ 
structed after calculations of the atom or its low-charged 
ions. The basis set contains inner core spinors, marked 
below by index “G” (which are occupied by the inert 
electrons and excluded from molecular calculations with 
GRECP as inactive, completely occupied states); va¬ 
lence, outer core (occupied by explicitly treated core elec¬ 
trons) and low-lying virtual spinors; all together they are 
marked by index “W”. The rest spinors, corresponding to 
high-energy virtual states, are marked by index “R”. The 
core states are only negligibly changed in the low-energy 
process under consideration (formation of chemical bond, 
low-energy excitation of atom, etc.). The completeness 
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condition for the {fnijm} basis can be formally written 
as 

1 = X! = Pc + Pw + Pr, (21) 

nljm 

Pc = ‘fnljmfhijmi (^2) 

nljm^C 


Pw — ^nljm^Pnijm^ (^3) 

nljm^W 

^R ~ V’n-ljmfnijm’’ (^ 4 ) 

nljm^R 

where Pc is the projector on the inner-core spinors, Pw is 
the projector on the outer-core, valence and low-energy 
virtual spinors, and Pr is the projection on the other 
(high-energy) states. For low-energy processes, which in¬ 
clude chemical bonding, low-lying excitations and those 
induced by weak external fields, one can usually neglect 
the high-energy states to study the properties of interest: 


P = {Pc + Pw + Pr)p{Pc + Pw + Pr) 

~ {Pc + Pw)p{Pc + Pw) (25) 
~ PcpPc + PwpPw = P^ p P^ , 


where p^ = PcpPc and p'^ = PwpPw- Here we have 
taken into account that the inner-core electrons need not 
be usually correlated to preserve high accuracy for the 
core properties in general, in contrast to the W-states. 
Therefore, the off-diagonal blocks PwpPc and PcpPw 
can be mostly neglected m- 

Due to Eq. (25) the expression (18) for (H) reduces to 
the following two terms: 

(H) « 


(^) 


w 


(26) 


where 
(A)^ = 

X) Pnljm,n'l'j'm' J Pn'l'j'm'^Pnljmdr , 

nljm^C,n' I' j'm' 


{A)^ = 

Pnlj7n,n'l'j^m' J Pn'l'j'm'' 

nlj'm^W,n'V j'm' ^W 

(27) 

For such properties as hyperfine structure constant, 
etc. considered here « 0 , i.e., direct contribution 

from the closed-shell core electrons can be ignored (and 
only their spin-polarization by open valence shells can be 
not negligible, see [SS]). In other cases {A)^ can be ob¬ 
tained from atomic calculation if one takes into account 
that the heavy-atom inner-core electrons are inactive in 
low-energy processes. Thus, for our case we have: 

(H) « {A)^ . (28) 

It is well known for heavy atoms and their com¬ 
pounds [HU [721 EHl EH] that the valence one-electron 
wave-functions and low-lying virtual states are propor¬ 
tional to each other in the vicinity of a nucleus. This 



FIG. 1. Large components of the 5 si/ 2 , 6 si /2 and 7 si /2 
spinors of Th for the configuration. The large com¬ 

ponents of 5 si/ 2 , 6 si /2 and 7 si /2 spinors in the core region 
are given in subfigure, where the scaling factor is chosen in 
such a way that the amplitudes of large components of these 
spinors are equal at Rc = 0.25 a.u. 


is due to overwhelming contribution of highly charged 
(even being shielded) heavy-nucleus potential as com¬ 
pared to other potentials from molecular environment, 
inter-electron interaction, etc. In figure one can see 
the large components of 5 si/ 2 , 6 si /2 and 7 si /2 of Th 
atom taken from a self-consistent field calculation of the 
Is^lp^Qd} configuration. Note that core radius Rc be¬ 
longs to the region of proportionality for the considered 
here core properties. 

One can introduce some reference functions for each 
combination of l,j: 


'Pijm{'r) = 


v!j{r)Sijm I 
Vfj{r)02j-i,jm J 


(29) 


The functions r](j^{r) will be determined such that they 
are equal to valence functions fij{r), gij{r) of a given 
atomic four-component spinor with the same Ij for r<Rc 
and equal to zero outside the (given core) region. Due to 
the proportionality property (see figure]^ it is not prac¬ 
tically important which of the W-functions, pYijjn{'^j is 
chosen for a given Ij or from which configuration it is 
chosen EZj. For example, for the case of Z = 0, j = 1/2 
of Th one can consider 5 si/ 2 , or 6 si/ 2 , or 7 si /2 as follows 
from Fig. For all functions from the W-diversity (i.e. 
for all nGW) 

‘pnljm{P) ~ ^nljm'Pljm{P) ; ^ ^ Pc : (30) 

where knijm are the proportionality (scaling) factors. 

Using Eqs. (20) and (30) we can rewrite Eq. (28) in 
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the following form: 


(A) « ^ Pnljm,n'l'j'm' j ^\iiij'rn'-^Vnljmdr 

~ 'Ylt Pnljm^n'Vj’m'^nljm^n’Vj'm' J j'm' 

nljm\n'l'j'm'^W 

ljm-,l' j'm' 

(31) 

where 


A 


ljm,l' j'm' — 


E 


nljm\n' I' j'm' 


Pnljm^n'l' j'm' knljm^n'l' j'm' ■ 

(32) 

Here is the W-reduced density matrix, in 

which the terms are summed up on the chosen principal 
quantum numbers n & W in contrast to a conventional 
one-electron density matrix from Eq. (fll 


The last expression in Eq. (31) means that for cal 


culation of a core property A it is sufficient to know 
some W-reduced density matrix as well as ma¬ 

trix elements of the operator over the reference functions 
{Tdijmii^}- One can interpret \\Aijm,i>j'm'\\ as a density 
matrix of an effective AlC state. The diagonal elements 
of the matrix are occupancies of the reference functions 
[78| . The nondiagonal elements between different Ij are 
“overlap occupancies (populations)” The latter can oc¬ 
cur in consideration of a molecule/crystal due to polar¬ 
ization of atomic orbitals in molecular/crystalline envi¬ 
ronment or in an atom placed in some external field, i.e. 
they reflect non-spherical distribution of electron density 
in the vicinity of nucleus of atom under consideration. 
We should stress that the diagonal and overlap popula¬ 
tions have meaning of “observable quantities” (though, 
in practice, their combination can be rather observed ex¬ 
perimentally). It means that different parts of the W- 
reduced density matrix can be obtained from different ex¬ 
periments (or their combinations). In particular, the di¬ 
agonal matrix elements can be extracted from the X-ray 
emission chemical-shift experiments. The non-diagonal 
matrix elements predetermine the value of effective elec¬ 
tric field which, in turn, can formally be extracted from 
the electron EDM experiments if we know the eEDM 
value. Finally, one can say that the AIC effective state 
in some sense is a more general term than the classical 
effective state term since the AIC effective state can in¬ 
clude overlap populations between different harmonics as 
is discussed above. 

The AIC concept described here can be applied both 
in the direct four-component calculation and in the two- 
step study, in which the four-component density matrix 
is obtained at the second stage of the procedure dis¬ 
cussed in the previous section. Actually, computation 
of the effective W-reduced density matrix is a special 
case of the recovery procedure when the equivalent ba¬ 
sis sets are constructed only from the reference functions 
{diijmir)} and some modification of one-center restora¬ 
tion is applied. In the present paper we report imple¬ 


mentation of the procedure and its application to calcu¬ 
lation of the W-reduced density matrix for the ThO and 
ThF'^' molecules and a number of core-properties: hy- 
p^rfine magnetic dipole constant (101, effective electric 
neld (i). the molecular-structure parameters of T,P-odd 
pseudoscalar-scalar electron-nucleus interaction (|^ and 
T,P-odd interaction of the nuclear magnetic quadrupole 
moment with electrons (|^. The code is interfaced to the 
DIRAC 12 [M] and MRCC codes. 


COMPUTATIONAL DETAILS 

To evaluate Egg, Wt,p, Wm and A|| in ThF+ we have 
applied the two-step method described above. The com¬ 
putational scheme used in the present paper is similar to 
that employed in nziiin] for calculation of ThO, where we 
have described and analysed the scheme in details (possi¬ 
ble sources of errors, importance of correlation treatment, 
importance of multireference approaches, applicability 
and convergence of multireference configuration interac¬ 
tion approaches, etc.). In all the calculations the Is—4/ 
inner-core electrons of Th were excluded from molecu¬ 
lar correlation calculations using the valence (semi-local) 
version of the GRECP [63] method. The main calculation 
was performed within the 38-electron two-component 
single-reference coupled-cluster method with single, dou¬ 
ble and perturbative triple cluster amplitudes, 38e-2c- 
CCSD(T). The calculation was perform using MBas basis 
set, generated in [21] with added h and i type functions, 
i.e., we used the (30,20,10,11,4,6,5)/[30,8,10,4,4,2,1] ba¬ 
sis set. For F we have applied the aug-ccpVQZ basis set 
[71] with two removed 5-type basis functions, i.e., the 
(13,7,4,3)/[6,5,4,3] basis set was used. To consider high- 
order correlation effects we calculated correlation correc¬ 
tion. For this we have frozen 20 outer core electrons 
(5s^5p®5(i^° shells of Th and Is^ shell of F) and per¬ 
formed two-component calculations within the coupled- 
cluster method with single, double, triple and pertur¬ 
bative quadruple cluster amplitudes, CCSDT(Q), and 
within the CCSD(T) method. We utilized the CBasSO 
atomic natural basis set which was generated using the 
same procedure that was used and described in [211 HI] 
and can be written as (35,29,15,10,7)/[6,8,5,3,2] for Th, 
and (13,7)/[4,3] for fluorine. The correction was calcu¬ 
lated as a difference between the calculated parameters 
within the CCSDT(Q) and CCSD(T) methods. In addi¬ 
tion, the basis set enlargement corrections to the con¬ 
sidered parameters were also calculated. For this we 
have performed: (i) scalar-relativistic CCSD(T) calcu¬ 
lation using the same basis set as used for the main two- 
component calculation; (ii) scalar-relativistic CCSD(T) 
calculation utilizing the extended basis set on Th (Lbas 
basis set (37,29,15,14,I0,10,5)/[22,17,15,14,10,10,5] gen¬ 
erated in [21] ). Corrections were estimated as differences 
between the values of the corresponding parameters. Fi- 
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nally, we have calculated the vibrational contribution to 
the considered core properties and molecule-frame dipole 
moment corresponding to zero vibrational level of the 
^Ai electronic state as a difference between the value 
averaged over zero vibration wave function and the non- 
averaged value at the given internuclear distance (3.75 
a.u., see below). The potential energy curve was calcu¬ 
lated at the 38-electron one-component CCSD(T) level 
with the LBas basis set. 


RESULTS AND DISCUSSIONS 

According to the 38-electron two-component CCSD(T) 
calculations the equilibrium internuclear distance in the 
^Ai state of ThF+ is 3.75 a.u. which agrees well with the 
experimental datum [35] , see table |T] In calculations of 
the parameters under consideration we have set R(Th“ 
F) to 3.75 a.u. 

Table nil lists the calculated values of effective elec¬ 
tric field along with the parameter of the T,P-odd 
pseudoscalar—scalar electron—nucleus neutral currents 
interaction, hyperfine structure constant, Wm parame¬ 
ter and g-factor for the ^Ai state of ThF+. It follows 
from table |ll| that the calculated value of EeS is sta¬ 
ble with respect to the electron correlation improvement 
and basis set enlargement. Similar to Ref. [29] using the 
size-extensive coupled-cluster calculations we have found 
that the outer-core electrons of Th contribute about 3.5 
GV/cm to EeS (similar value was found in Ref. [55] for 
ThO) and -161 ^^-MIIz to A||. Thus, if one performs 
18-electron rather than 38-electron calculation the outer- 
core contributions should be taken in mind. According 
to our calculations the spin-orbit contribution from the 
valence electrons to Aeff (about 1 GV/cm) is almost neg¬ 
ligible in the case of ThF+, in contrast to ThO where 
it is about 10 GV/cm [SS])- According to table [IT] in 
view of the extensive analysis of uncertainties performed 
in Ref. [55] we suggest that the theoretical uncertainties 
of Fleff, Wt,p and A|| are within 7%. Unfortunately, A|| 
(Eq.fiol) is yet unknown experimentally for ^^®ThF+ and 
it cannot be used currently to check the value of Ef.s and 
other considered properties. However, we have shown in 
Ref. [53] for the ground state of PbF molecule that the 
used computational scheme is rather accurate: the cal¬ 
culated value of A||(PbF) agrees with the experimental 


TABLE I. Equilibrium internuclear distance Re, harmonic vi¬ 
brational wavenumber cUe and vibrational anharmonicity uj^Xe. 
for the ^Ai state of ThF^. 


Method 

Re, a.u. 

LOe, cm ^ LOeXe, Cm ^ 

MRci-tQ/so, gg 

3.76 

655.6 — 

CCSD(T), this work 

3.75 

658.4 1.9 

Experiment, [43] 

3.74(4) 

658.3(10) — 


datum [ST] within 2%. Finally, it should be noted that 
the estimation made in Ref. [26] for E^s (90 GV/cm) is 
more than twice overestimated (similarly strong overes¬ 
timation was also found for E^g in PtH+, see [55]). In 
a like manner the estimations for the Wt,p and Wm pa¬ 
rameters made in Refs. [55] [55] and based on EeS from 
Ref. [55] are also about twice overestimated. 

The effective electric field in the ^Ai state of ThF"*" is 
about two times smaller than the E^g in the ^Ai state of 
ThO (81.5 GV/cm, see [55]) because of a smaller mixing 
of s and p orbitals. We can give the following explana¬ 
tion. In the naive ionic model ThO can be considered 
as Th+^ and 0“^, ThF+ can be considered as Th+^ and 
F“^. This agrees with the fact that the molecule-frame 
dipole moment of ThO is about 1.5 times larger than the 
dipole moment of ThF"*" with respect to the Th nucleus 
(see Table II). This leads to higher effective negative elec¬ 
tric charge on oxygen in ThO than on fluorine in ThF+. 
Both ThO and ThF+ have two unpaired electrons. They 
are non-bonding and localized on Th so that Th has 
configuration in both cases, where a is mainly the 7s 
atomic orbital of Th and 6 in mainly 6d atomic orbital of 
Th. The unpaired electrons of Th feel a stronger electric 
field in ThO than in ThF+. This leads to higher polar¬ 
ization of the unpaired electrons in the case of ThO, i.e., 
stronger s—p mixing of a-state {S state has no practical 
interest for E^g here due to far smaller amplitude of 6d 
in the core region than 7s). The leading contribution to 
Egg is roughly proportional to 


Cr,Crp{7s\Hd/dg\7p) , 

where Eq. 0 is used. The matrix element is mainly 
accumulated near the Th nucleus and l);C 7 p are 

the corresponding MO LCAO coefficients of the atomic 
Th orbitals in the hybridized molecular one. As a con¬ 
sequence the effective electric field should be expected 
notably larger in ThO vs. ThF+. On the other hand the 
smaller polarization of open-shell cr-state leads to higher 
s—character of the orbital and the hyperfine structure 
constant in ^ Ai of ThF+ is bigger than in ThO (see table 
E- Note, that the hyperfine structure constant behaves 
“inconsistently” with respect to the effective electric field 
in the present case. 

We have applied the AIC theory described in the pre¬ 
vious section to the case of ThF+ and ThO. To set the ra¬ 
dial reference functions given in (29) we have used 

7s, 7p and 6fi functions from calculation of the 7s^7p^6d^ 
configuration of Th and have evaluated the W-reduced 
density matrix defined by Eq. ( [3^ from the 

molecular density matrices obtained within the CGSD 
approach. Note that {T-Lijm} (29) coincide with the 7s, 
7p and 6d functions within some radius Rc (here we set 
Rc = 0.25 a.u.) and are zero outside the radius. For 
brevity we will designate the selected reference functions 
as 7 si/2,i/2 : etc. The operator of hyperfine interaction 
mixes the states with same parity and to; the diagonal 









TABLE 11. The calculated values of the molecule-frame dipole moment (d), effective electric field (iSeff), parameter of the 
T,P-odd pseudoscalar—scalar electron—nucleus neutral currents interaction {Wt,p), parameter of T,P-odd MQM interaction 
(Wm), hyperfine structure constant (A||) and g-factor (G||) of the ®Ai state of ThF^ compared to the corresponding values of 
ThO from Ref. |28i I29 | using the coupled-cluster methods. 


Method 

d *, Feff, Wt,p, 

Debye GV/cm kHz 

Wm 

lO^^Hz 
e cm^ 

A||, 

^•MHz 

mn 

G|| 

38e-2c-CCSD 

2.69 

35.5 

48 

0.87 

-4214 

0.039 

38e-2c-CCSD(T) 

2.66 

38.1 

51 

0.90 

-4164 

0.033 

correlation correction 

0.07 

0.0 

0 

-0.01 

13 

0.001 

basis set correction 

-0.01 

-0.6 

-1 

-0.02 

-14 

— 

vibr. contribution 

0.03 

-0.1 

0 


2 

— 

FINAL(ThF+) 

2.74 

37.3 

50 

0.88 

-4163 

0.034 

FINAL(ThO) 

4.23 

81.5 

112 

1.66 

-2949 0.007 


(see Ref. [28l[29]) _ 

* The dipole moment is calculated with respect to Th nucleus. 
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matrix element of the operator for \ljm) state is opposite 
by sign to the diagonal matrix element for \lj — m) state. 
Thus, the diagonal contribution to the mean value of the 
operator is defined by the difference between the density 
matrix elements, Aijmjjm - and the diag¬ 

onal hyperfine operator terms for \ljm) states. For the 
most important elements of W-reduced density matrix 
for the ThF+ ^Ai state we have: 


Ay y - Ay y 

7Si/2,1/2.7si/2,1/2 7si/2, _ 1/2,7Si/2, _ 1/2 

A^ ^ - A^ ly- 

‘Pl/2,l/2)'Pl/2,l/2 ‘Pl/2,-l/2’'Pl/2,-l/2 

6ti3/2,3/256^3/2,3/2 6^3/2,-3/2 56^3/2,-3/2 


“h ^'7' 

‘ ^l/2,l/2'>‘Pl/2,l/2 ‘^1/2,-1/25'Pi/2,-1/2 

-t" ' 7 ' 

7Pl/2,l/2)7si/2,l/2 7pi/2,_i/2 57Si/2,-l/2 


-0.99 

-0.47 

0.88 

0.105 
0.105 . 


(33) 

For the most important elements of W-reduced density 
matrix for the ThO ^Ai state we have: 


Ay y - Ay y = -0.72 

'®l/2.1/2)‘Sl/2,l/2 <^1/2,-1/2 5'''51/2,-1/2 

Aiy- ^ - Aiy- = -0.37 

'Pl/2,l/2> 'Pl/2,1/2 ‘Pi/2,-1/2 I 'Pi/2,-1/2 

A^ ^ — A^ ^ = 0.67 

603 / 2 , 3 / 25603 / 2 , 3/2 603 / 2 ,- 3 / 25603 / 2 ,- 3/2 


7Sl/2,l/2,7pi/2,i/2 


7p 


1 / 2 , 1/2 


,7si 


/2,l/2 


A.- =0.238 

7Si/2,-l/257Pi/2,-l/2 

7 % =0.238 

'Pl/2,-1/25'^1/2,-1/2 


(34) 


Thus, in terms of the reference functions the effective 

configuration of unpaired electrons of Th in ThF+ is 
— 1 . 0 —^ 0 . 5 —^ 0.9 — 0 . 7 —^ 0 . 4 —^ 0.7 

7s 7p 6d , while in ThO it is 7s 7p 6d . The 

leading matrix element of the HFS operator in the basis 
of reference functions is between the 7s functions. The 
ratio of effective occupancies of 7s « 1.4. This explains 
the ratio of the HFS constants given in table |llj 

The operator of effective electric field ([^ mixes states 
of opposite parity with the same m; matrix element of 
the operator between \ljm) and \l'jm) has the same sign 
as the matrix element between \lj — m) and \l'j — m) 
state. Thus, the mean value of the operator is defined by 
combination of the types Aijmpjm + j-m of the 

W-reduced density matrix \ and matrix ele¬ 

ments of ifeff operator between \ljm) and \l'jm) states. 
Most important of the combinations of ||Aij„i^;/j/m'|| 
matrix elements for ThF+ and ThO ^Ai are given in 
Eqs. (33 34). From the equations one can see that the 
W-reduced overlap population between the 7s and 7p 
functions in ThO is twice larger than that in ThF+. This 
explains the appropriately larger Eeff in ThO, see table 

im 


Note that the W-reduced density matrix \ \Aijm,i'j'm'\\ 
(and minimal number of core-property matrix elements 
over the reference W-reduced functions) can be con¬ 
sidered as a pretty concise description of the effective 


atoms in compounds state which is appropriate for “al¬ 
most quantitative” calculation of the mean values of core¬ 
properties under consideration. 


CONCLUSION 

The parameters Egg, Wt,p and Wm which are re¬ 
quired to interpret experimental measurements on the 
^Ai state of ThF+ in terms of fundamental quantities 
are calculated. The value of E^g in ThF+ was found to 
be notably smaller than that in ThO. A quantitative ex¬ 
planation is given. On the other hand, i?eff(ThF+) is 1.6 
times bigger than the effective electric field in the HfF+ 
cation da [ig which is under preparation for electron 
electric dipole moment search [T5] . 

In the present paper we have implemented the concept 
of atoms in compounds and applied it to calculate the 
W-reduced density matrix for description of the effec¬ 
tive state of Th in ThF+ and ThO. This matrix contains 
“sufficient information” to evaluate such physically ob¬ 
servable properties as hyperfine structure constant, etc. 
whereas the conventional density matrix is excessive here. 

According to our preliminary study, the electronic 
spectrum of ThF'*" is more dense compared to the ThO 
molecule. Its accurate theoretical investigation requires 
inclusion of quadruple cluster amplitudes as shown in 
Ref. [ig and is also found in our preliminary study of 
^Ai-^S transition energy. We plan to investigate it in 
our future study of ThF+ elsewhere. 
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